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Abstract. The purpose of this paper is to obtain exact solutions 
of the Einstein field equations describing traversable wormholes 
supported by phantom energy. Their relationship to exact solu- 
tions in the literature is also discussed, as well as the conditions 
required to determine such solutions. 
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1. Introduction 

Wormholes may be defined as handles or tunnels in the spacetime 
topology linking different universes or widely separated regions of our 
own universe |1 . Renewed interest in the subject is due in part to 
the discovery that our universe is undergoing an accelerated expansion 
12 E], that is, a > in the Friedmann equation a/ a = — ^(p + 3p). 
(Our units are taken to be those in which G = c = 1.) The accel- 
eration is caused by a negative pressure dark energy with equation of 
state p = —Kp, K > ^ and p > 0. A value of > | is required for 
accelerated expansion; K = 1 corresponds to a cosmological constant 
Of particular interest is the case K > 1, referred to as phantom 
energy. For this case, p + p < 0, in violation of the null energy condi- 
tion. (The null energy condition requires the stress-energy tensor Ta/s 
to obey Tapk'^k^ > for all null vectors.) It is well known that the 
violation of the null energy condition is a necessary condition for the 
existence of wormholes p] . In this context such matter is usually called 
exotic. Phantom energy could therefore automatically qualify 
didate for exotic matter, except for one problem: the notion of dark 
or phantom energy applies to a homogeneous distribution of matter 
in the universe, while wormhole spacetimes are necessarily inhomoge- 
neous. Fortunately, the extension to spherically symmetric inhomoge- 
neous spacetimes has been carried out. (See Ref. j3] for details.) 
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Returning to the definition of phantom energy, we saw that the con- 
dition K > 1 results in a violation of the null energy condition. So if 
K < 1, the null energy condition is met and turns out to have a direct 
effect on the so-called fiare-out condition. We shall return to this point 
in Section m 

A recurring problem in the general theory of relativity is finding exact 
solutions to the Einstein field equations. In the case of phantom-energy 
supported wormholes several solutions already exist in the literature 
El El- The strategy in this paper is to start with a general line 
element, together with the above equation of state, and to determine 
the conditions required to obtain explicit exact solutions, both old and 
new. Two new solutions are discussed. 

One problem is that not everyone uses the terms exact and explicit 
in the same sense. Thus an "exact" solution may contain an intractable 
integral, while an "explicit" solution may contain functions defined only 
implicitly. Our primarily interest is therefore centered on elementary 
functions, rather than arbitrary functions. For present purposes these 
may be defined as functions of a single variable built up by using that 
variable and constants together with a finite number of algebraic opera- 
tions, composition, forming trigonometric functions and their inverses, 
and constructing exponents and logarithms. 

While the derivative of an elementary function is elementary, the in- 
tegral may not be. For example, J e^ dx is not an elementary function, 
or, as it is often expressed, the integral cannot be written explicitly (or 
in closed form or in finite terms.) General criteria for integration in 
finite terms can be found in Refs. jHUll- 

We have similar requirements for the solution of differential equa- 
tions. We are interested in finding solutions that can be expressed ex- 
plicitly in terms of elementary functions, as opposed to infinite-series 
or numerical solutions. 

The solutions in Refs. [SIEIIZI mentioned above are examples of exact 
solutions in the sense defined here. 

2. The problem 
Consider the general line element 
(1) ds^ = -e^^Wdt^ + e^^Wrfr^ + r\de^ + sin^^ dc/.^), 

where $ and a are functions of the radial coordinate r. The function 
$ is called the redshift function. We require that e^**^^-* never be zero 
to avoid an event horizon. The function a has a vertical asymptote at 
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the throat r = tq: 

hm a{r) = +00. 
The reason is its relationship to the shape function b{r): 

(2) e^"^ = ^TTT- 

X _ Mr) 

r 

It follows that 

(3) 6(r) = r(l-e-2"M). 

The shape function determines the spatial shape of the wormhole when 
viewed, for example, in an embedding diagram. To obtain a traversable 
wormhole, the shape function must obey the usual fiare-out conditions 
at the throat P]: b{ro) = ro,6'(ro) < 1, and b{r) < r. Another require- 
ment is asymptotic flatness, that is, h{r)/r — > as r ^ cxd. Because of 
the spherical symmetry, the nonzero components of the stress-energy 
tensor are Too = p{.f)^ Tn = p{r), and T22 = T33 = pt{r), the transverse 
pressure. The components of the Einstein tensor in the orthonormal 
frame are given next [TUj : 

(4) Gtf = -e-^^Wa'fr) + ^(1 - e^^^M), 

(5) Gff = -e-2°M$'(r) - ^-(1 - e'^^M), 



(6) Gqq - G^^ 

= e-2"(") ( $"(r) - <l>'(r)a'(r) + [^'(r)]^ + 



r) Ci'i^] 



From the Einstein field equations G^^ = SttT^^ and the equation of 
state p = —Kp, we have G^^ = Svrp and Grf = 87r(— fCp), giving us the 
following system of equations: 

2 . 1_ 

„2 ' 



(7) Gii = SvrT,, = Svrp = -e-2"Ma'(r) + -(1 - e'^^M^ 



(8) Gff = SnTff = 8n{-Kp) = -e~'"('")<l''(r) - ^(1 - e-'"^")). 



Substitution yields 
-e-2"Ma'(r) + ^(l-e-2°«^ 
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After rearranging the terms, 

(9) Ka\r) = -$'(r) - ^(e'"^'') - 1){K - 1). 

This equation shows the close relationship between $'(r) and a'{r) 
and hence between $(r) and a{r). Since a{r) — > +oo as r — > ro+, there 
is a distressing tendency for e^*^''^ to go to zero as r — > ro+. So, while 
the existence of exotic matter does help to satisfy a basic requirement, 
the equation of state makes it very difficult to obtain an exact solution 
without an event horizon. 



3. The redshift function 

Returning to Eq. one way to solve this equation is to insert the 
redshift function "by hand." One possibility is $'(r) = 0, resulting in 
$ = constant. This equation is readily solved and leads to 

(10) e^-M = ^-ttTk- 

This is the solution in Ref. [Hj. 

The only other possibility is $(r) = |ln(ri/r), for some constant ri, 
which allows the solution of Eq. Q by separation of variables, that is, 
by factoring 1/r. This approach yields 

(11) e^^M 



This is the solution in Ref. |7j. 

It follows that to obtain an exact solution, $'(r) must depend di- 
rectly on a{r) and a'{r) and therefore indirectly on the shape function. 
This dependence may be expressed as $'(r) = F[a{r)]a'{r), for some 
elementary function F. 

4. The general case 
Suppose we write Eq. © as follows: 

(12) -Ka'(r) - — (e^^M - 1)(K - 1) = $'(r). 

2r 

By the above discussion, $'(r) must have the form 

(13) $'(r) = F[a{r)]a'{r). 
So by Eq. 

Ka'ir) 1 F[a{r)]a'{r) 

(14) -3-73 — --— (A-i;- 



Q2a{r) _i 2r e^"^'') - 1 
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To obtain an exact solution, we must be able to solve Eq. (|T^ in a 
closed form and to express the integral of $'(r) in finite terms. Finding 
$(r) in an exact form does not, of course, guarantee the absence of an 
event horizon. For example, if F = — 1 in Eq. ()13|) . then $(r) = 
—a{r) + C and e^**^^'-* = at the throat, which is a typical outcome! 
The simplest way to meet all these requirements is by letting 

^[«(^)] = 

Then, Eq. (fTlj) becomes 

Ka'(r) 1,^^ , Ka'(r) 
-{K - 1, 



with initial condition a{ro) = +oo. The solution is 



Integrating —Ka'{r)e~'^°'^'^^ and substituting the expression for e^"*-'"-' 
yields 

1 / - ^ 
$(r) = - In C - 

C is the constant of integration, which needs to be determined from 
the junction conditions. So, the line element is 

(15) ds^ = -C(-] de + -———-^dr^ + r\de^ + siYi^ed(P^). 
\roJ In(-)v^-i)/^ 

This solution has most of the required features. For example, since 
b{r) = r(l - e-2"('-)) by Eq. ©, we have 



6(r) 



1 - In 



K-l 
K 



It is easily checked that &(ro) = tq, Vijo) = 1/K < 1, and 6(r) < r. 
(Observe that if < 1, then &'(ro) > 1 and the flare-out condition is 
no longer satisfied.) 

Unfortunately, the resulting spacetime is not asymptotically flat since 
b{r) eventually decreases. Accordingly, the wormhole material must be 
cut off at some r = a and joined to an external Schwarzschild space- 
time. A natural choice for r = a is the value for which b{r) becomes a 
maximum. From the critical value of b'{r) = 0, we get a = roe^^^^~^\ 
At first glance this does not look like a large distance. According to 
Ref. [TT], however, K is likely to be very close to unity. 
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Matching our interior solution to the exterior Schwarzschild solution 



(16) ds'^ 



1 - 



2M 



2M 



dr^ + r^{de^ + sm'^9d^'^) 



at some r = a requires continuity of the metric. As noted in Ref 
since the components g^g and g^^ are already continuous due to the 
spherical symmetry, one needs to impose continuity only on the re- 
maining components at r = a: 



and 



fi'rr(int) (o) — fi'rr(ext) {o) 



9tt{mt) l^J — 9it{cxt) [(^) 

for the interior and exterior components, respectively. These require- 
ments, in turn, imply that 



and 



In particular. 



bmt{a) = fecxt(a)- 

1 



1 _ ^ 1 _ 2M • 



So we need to determine M = \b{a), the total mass of the wormhole 



for r < a: 



M = -h{a) 



1 - In 



Choosing a = tqc 



1/{K-1) 



we get 



M = -roe^-i ( 1 



Returning to $(r), we now have 



or 



To 



c 



K-1 



1 

K 



2M 



1 - 2M/a 



(a/ro)^'-i- 

For a = rQe^^^^~^\ C = l/{Ke). So the line element becomes 
(17) 



ds^ = - 



Ke V ro 



dt^ + 



dr^ + r\d9^ + sm'^ 9, 



While the metric is continuous on the junction surface r = a, the 
derivatives may not be. This behavior needs to be taken into account 
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when discussing the surface stresses. The following forms, proposed by 
Lobo 0112], are suitable for this purpose: 



Alia 



2M / 6(a) 



and 




iira 



It is clear that the surface stress-energy a is zero. The surface tangential 
pressure V turns out be positive. To see this, consider $'1, the "jump" 
at r = a = roe^^^^~^^: lim^^a+P is positive since we are dealing with 
a Schwarzschild spacetime ^]. Substituting $'_(a) = limj,_^a_ $'(r) 
shows that lim^^^.P = 0. 

An important consideration affecting the traversability is the proper 
distance £{r) from the throat to a point away from the throat: 



£(r) 



dr 



which is finite; in fact, i{ro) = 0. Unless K is extremely close to unity, 
i{r) is not going to be excessively large. For example, if i^' = 1.1 and 
r = 2ro, then i{r) ^ 7.1ro. 

A final consideration is the time dilation near the throat. Let v = 
di/dr, so that dr = di/v (assuming that 7 = a/1 — [y/cY 1.) Since 
dt = e^'^'^^dr and dr = e^^^^dt, we have for any coordinate interval At: 

rtb rh Ap rn 1 

At= / dt= I e-*W-= / ie-*We"Mdr. 

Jta 



Going from the throat to r, we get 
At = 



Ke I frQ\^-'^ dr 
r 



^ln(^)(^-i)/^ 

yKe dr 



< 



which is also well behaved near the throat. 
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5. Other solutions 

Returning to Eq. (fTT^ . the choices for F appear to be severely hmited. 
At least one other possibility is F[a{r)] = -2K / {e^''^''^ + 1). So by 
Eq. (ini), 

e2a(r) _i 2r e'^"('') - 1 ' 

This equation can also be readily solved to yield 
(18) 

/ N 2{K-1) , 

ds^ = -C[-] de + ^^..^ ^^.^ dr'' + r\de'' + sin^^ dcf?). 

This solution has some of the same features as the solution in the 
previous section: h{r) satisfies the flare-out conditions and attains a 
maximum value at 



a = ro 



K 

2K \~ 



It is readily shown that if r = a is the junction surface, then 



M = tq- 



K -I f 2K \^ 



2K -1 \ 2K -1 
and 

(2K-ir^ 

{2KY^ ■ 

The cases considered so far may very well exhaust the possibilities for 
getting exact solutions without event horizons. For example, F[a(r)] = 
-l/Ve2"M - 1 in Eq. (O leads to 

^ r a\r)dr ^ _tan- Ve^-M - 1 + C, 
J Ve2°(^') - 1 

which is finite at r = tq. So there is no event horizon. Unfortunately, 
the differential equation (fT^ leads only to an implicit solution for e"^''-': 

_ ^ [ln(e2"(") - 1) - 2a(r)] - ^^^^Inr 



- (-2 tan-i Ve^^C'-) - 1 - , ^ ) + C. 



This equation cannot be solved explicitly for e°'^'^\ 

The choice F[a(r)] = l/(e^°'^'^^ — 1)^ is even worse. Not only do we 
get a strictly implicit solution, but the resulting metric has an event 
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horizon. Other plausible choices, such as 

= , ^ ,1 — — for n>2 

or 

are even more complicated and do not yield explicit solutions. 

Specifying a{r) and solving for $(r) leads to similar difficulties; judg- 
ing from Eq. (jl2p . 

$'(r) = -Ka'ir) - — (e^"^ - 1)(K - 1), 
2r 

if 

Jro 2r 

is finite, then 

$(r) = -Ka(r) - [ — (e^^M - 1)(K - l)dr 

is likely to lead to an event horizon at r = rg. But, as we have seen, 
$(r) can be determined in certain special cases, such as Eqs. (fTT1|l and 

m- 

Another special case is discussed in Ref. |j5j: if p(r) = pq, a constant 
for To < r < ri, then Eq. (|7j) becomes 

Svrpo = -e-2"Wa'(r) + hi - g-^^W). 

So 

8w' = 2re-2"('')a'(r) + (1 - e-'"^'')) = h'{r) 

by Eq. Q. This equation yields 6(r) and hence a{r) and $(r), although 
the details are complicated. 



6. CONCLUSION 

It is shown in this paper that for wormholes supported by phantom 
energy the only specific choices for the redshift function $ that lead to 
explicit exact solutions without an event horizon are $(r) = constant 
and <I>(r) = iln(ri/r). Otherwise $(r) must depend on a{r) in such a 
way that $'(r) = F[a(r)]a'(r), for some elementary function F. The 
choices for F are severely limited. Two new solutions are obtained. 
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